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I review recent developments in the study of strongly interacting field theories with non-zero chemical potential 
fi. In particular I focus on (a) the determination of the QCD critical endpoint in the (/i, T) plane; (6) superfluid 
condensates in Two Color QCD; and (c) Fermi surface effects in the NJL model. Some remarks are made 
concerning the relation of superconductivity with the sign problem. 



quark matter is unstable with respect to diquark 
condensation (qq) ^ 0, resulting in a ground state 
with color superconducting properties |^ . Model 
calculations suggest that the BCS gap A at the 
Fermi surface may be as large as 100 MeV [|, 
comparable with the constituent quark scale, im- 
plying significant consequences for the physics of 
compact astrophysical objects Q|. In this case 
simulations of QCD are to date impracticable, for 
reasons I will review below; however, diquark con- 
densation and Fermi surface effects can be studied 
by lattice techniques in certain models. 

1. Why is /i 7^ so difficult? 

For a vectorlike gauge theory in Euclidean met- 
ric the introduction of a quark chemical potential 
breaks the 75-hermiticity of the Dirac operator: 
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Figure 1. Schematic view of the QCD phase dia- 
gram. 



Fig. |i| summarises our current knowledge of 
the QCD phase diagram in the plane of temper- 
ature T and quark chemical potential /i. Last 
year Shinji Ejiri j|] reviewed QCD simulations at 
T 7^ 0, ie. along the vertical axis. In my talk I 
wish to discuss what can be done in the interior of 
the plane. Significant progress has been made in 
the region to the upper left of Fig. |l] where both fj, 
and T differ from zero, which is also of direct phe- 
nomenological interest for heavy-ion collisions. I 
also wish to cover the lower right region describ- 
ing cold dense strongly-interacting matter. In re- 
cent years there has been intense theoretical ac- 
tivity in this region driven by the possibility that 
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This implies that eigenvalues of p are no longer 
pure imaginary and hence not related to each 
other by complex conjugation, in turn implying 



detMifi) ^ dctM*{fi) = detM{-fi). 



(2) 



Therefore the functional measure is no longer pos- 
itive definite. In principle the determinant can be 
factorised into a modulus p and a phase 0, and the 
phase included with the observable O in Monte 
Carlo simulations via 



(0) ^ (0e^^),/(e^^)p. 



(3) 



Unfortunately fairly general arguments suggest 
that (e*"^)p oc e^^, where V is the system vol- 
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ume. Acquiring sufRcient statistics therefore be- 
comes exponentially difficult as the thermody- 
namic limit is approached. This is known as the 
Sign Problem, and has plagued the study of 7^ 
since its inception. However, we should not be 
so surprised, since generic problems in the quan- 
tum theory of TV objects require A^! ^ differ- 
ent complexions or wavefunctions to be examined. 
Perhaps a more productive way of phrasing the 
question would be why is vacuum QCD so easy? 

Given this difficulty there are two routes for- 
ward. Firstly, one can perform a QCD simula- 
tion at /I = and attempt to analytically con- 
tinue the results to /i 7^ 0. This can be done 
either by calculating terms in a Taylor expansion 
about = 0, exemplified, for instance, by the 
calculation of baryon number susceptibility xb = 
y-i^^inZ/VU^o il, or by directly reweight- 
ing configurations, as in the 'Glasgow method' 
developed by Ian Barbour and collaborators ||]. 
In the former case the prospects are limited by a 
finite radius of convergence, dictated by the pres- 
ence of a critical point in the (/i, T) plane. In 
the latter case the effectiveness is limited by the 
requirement of maintaining a reasonable overlap 
between the trial and true ensembles, which once 
again becomes exponentially hard in the thermo- 
dynamic limit 10]. Both methods have their best 
chance of succeeding at T > 0; at T = the be- 
haviour is horribly non-analytic, since the ground 
state does not change as /i increases out to the 
onset of nuclear matter at fio — 307MeV. As we 
shall see, the overlap problem is also considerably 
more tractable at T > 0. 

The second approach is to throw up our 
hands and use a real measure detMM* = 
detM(/Lt)M(— /z). Physically this has the effect of 
introducing conjugate quarks q'^ which have pos- 
itive baryon charge but transform in the conju- 
gate representation of the gauge group [||, lead- 
ing to the possibility of gauge invariant qq'^ bound 
states. In real QCD this is a disaster since the 
lightest such state is degenerate with the pion, 
being in effect a 'Goldstone baryon'. As /i is 
raised the onset of nuclear matter would thus 
be expected at fiQ « rather than the con- 

stituent quark scale E « m^/i- There are mod- 
els, however, such as Two Color QCD ^ and 



QCD with non-zero isospin chemical potential 
/i/ cx /ill — A*d where qq Goldstone baryons 
are a feature, not a problem Another class 
amenable to this approach are NJL-like models, 
where qq'^ states don't couple to the Goldstone 
mode fl^ and hence remain at the constituent 
scale, so the use of a real measure doesn't dam- 
age the physics. The common feature of all these 
models is that the diquark condensate which po- 
tentially forms at large is gauge singlet, and 
hence describes a superfluid, as opposed to su- 
perconducting, ground state. 

2. Progress at T > 

Fodor and Katz ||l3| have made spectacular 
progress in reweighting at T > 0. The basic for- 
mula is 

Z[a] = j DUeyi]i{-Sbos[U\ao\)detM[U;ao\ x 

f . A o rrr IX detM[f7;a] 1 
r"P(-^^^-[^'"'""])detM[[/;ao]/ 

where the parameter set a — {(3, m, /i} but im- 
portance sampling is performed using a different 
set ao = a — Aa, with chosen either zero or pure 
imaginary to keep the determinant real. Usually 
reweighting is only effective if Aa is small enough 
to maintain a good overlap between trial and true 
ensembles. Fodor and Katz' insight is that along 
the crossover/coexistence line between hadronic 
and quark-gluon plasma phases the true ensem- 
ble cannot alter very much, since it must contain 
contributions from both phases; hence the overlap 
between ensembles along the line remains high. 
They therefore generalise the Glasgow method 
by reweighting with both A/i, A/3 / in order 
to stay on this line. In practice PdfJ-) is iden- 
tified via the real part of the lowest Lee- Yang 
zero Poy which approaches the real axis as V in- 
creases. Now, the nature of the transition be- 
tween hadronic and QGP phases can be deter- 
mined by the volume scaling of Im/3o: 

1- T o f 7^ 0; crossover, 
hm Im/3o <, r, n. . , \ (4) 
y^oo I = order transition. 

The results for 2 + 1 flavor QCD with mu.d^ = 
0.025, TTisO = 0.2 on volumes up to 8"^ x 4 are 
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Figure 2. Reweighting prediction for the quark- 
hadron transition (note /x_b = 3/x). 



charged particle yields and strangeness enhance- 
ment. QCD-TARO have looked at the response 
of hadron masses to a change of both baryon 
and isospin chemical potentials They find 

d^rriTT /dfi^ > 0, and with a much larger value 
in the QGP phase, indicative that the pion is 
no longer a Goldstone mode at high T; by con- 
trast 9^m^/9/i| < and is much larger in the 
hadronic phase, consistent with pion condensa- 
tion in cold dense isospin-asymmetric matter JlOf . 
Finally Ejiri has determined the curvature of the 
critical line d"^ jSdlJi) / d^^ by measuring the shift 
in the peaks for {ip^) and Polyakov loop suscep- 
tibities jl^. For quark mass ma = 0.2 on 16^ x 4 
he obtains a value ~ —1.5(4); assuming a criti- 
cal temperature Tc ~ 170MeV and a /3-function 
given by its one-loop value, this yields a value 
rc(Ai = 242MeV)« 150 MeV, not too far from (|). 



shown in Fig. ^ The physical scale is set us- 
ing separate simulations at /3c (0) and /3c (mb), 
which yield a critical temperature T(.{^. = 0) = 
172(3)MeV. The main prediction is for the criti- 
cal endpoint of the first order line, at 

= 160(4) MeV ; /i^ = 242(12) MeV. (5) 

The result for /ig is unexpectedly large, and 
might undermine proposals to observe it directly 
at RHIC iQ. However, the endpoint is expected 
to move to the left as the chiral limit is ap- 
proached, so one should consider (||) more as a 
proof of principle at this stage. What should be 
stressed is that currently the results are limited 
by cpu resources rather than overlap problems as 

— > oo; it is certainly possible, therefore, that 
lattice methods can ultimately obtain an accu- 
rate determination of {he,Te)- 

There have also been developments in Taylor 
expansion methods, which can deliver informa- 
tion on a wider range of observables. Gavai and 
Gupta have reexamined quark number suscepti- 
bilities in quenched QCD ||l^, in particular find- 
ing that Xs increases markedly across the transi- 
tion to QGP, and approaches Xu,d for T ^ 2Tc, 
though all three depart significantly from ideal- 
gas values. These quantities are testable in ion 
collisions via eg. event-by-event fluctuations in 



3. Two colors matter 

For gauge group SU(2), the fermion determi- 
nant is trivially real since all matter representa- 
tions are either real, or pseudoreal (in which case 
detM = detT2M*T2). Moreover, if we specialise 
to the case of TV staggered fermions in the funda- 
mental representation, the usual U(A^) (g) lJ{N)s 
global symmetry at m = /i = is enhanced, via 

X, = {Xe,-x'Ir2) ; X,= ( ^° ) (6) 

to X ^ VX, X ^ XV\ V e U(2iV). For ad- 
joint quarks the same relation holds with the T2 
factors replaced by 1. Now, at /j, = we expect 
this global symmetry to be spontaneously broken 
by a chiral condensate (xx) 0. For fundamental 
quarks the breaking pattern is U(2iV) — + 0{2N) 
yielding N{2N + 1) Goldstones; for adjoint the 
pattern is V(2N) -> Sp(2iV) with N{2N - 1) 
Goldstones [Q. This differs in detail from the 
corresponding patterns for continuum fermions 
||Tl| ]. The important point is that besides the 
usual mesonic qq states, in general some of the 
Goldstones must be qq baryons, as anticipated in 
Sec. |[ It is also possible to identify diquark con- 
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densates for both 2 and 3 representations: 



992 

qqs 



(7) 



where in the 3 case the antisymmetric e acts on 
flavor. These are related to (xx) by a U(27V) rota- 
tion and are gauge invariant. Since, however, (^) 
is not invariant under the original U(1)_b, baryon 
charge is no longer a good quantum number. 

It is possible to treat TCQCD analytically by 
ignoring all excitations except the Goldstones, 
and writing down an effective action in the spirit 
of chiral perturbation theory (xPT) , in which the 
physical parameters are (xx) ^^'^ ™7r sd, fi ^ 0. 
Remarkably, the ^-dependent terms in the chiral 
Lagrangian are completely determined in terms 
of these two parameters by the global symmetries 
. At leading order in xPT a second order on- 
set phase transition is predicted at the rescaled 
chemical potential x = — 1: 
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(8) 



where the baryon charge density ub — {ip^/oip) 
and the subscript denotes values at /i = 0. The 
high-/z phase is a superfluid which forms as a re- 
sult of a Bose- Einstein condensation of weakly in- 
teracting diquark bosons. 

Let us briefly discuss the measurement of {qq). 
In a flnite volume the technicalities are identical 
to those involved in measuring the chiral conden- 
sate. A diquark source term JX*'^'''2X is introduced 
and the action rewritten in a Gor'kov basis MM'- 



]T2 \M 



JT2 



c^ix.x'l 

whence 

Z[]J] = j P[/Pf(2^[C/,j,j-])e 
The condensate is then defined by 

{qq{j)) 



X" \ _ 

X 



= .4* (9) 



1 d\nZ 1 , 
V dj 2V^ ' 



(10) 



(11) 



Of course, since j is not a physical parameter one 
wants the j ^ limit, but any favourite method 



for probing {ip'ip{'m = 0)) can be used, such as 
direct inversion of „4(j) followed by extrapolating 
j using a Banks-Casher relation on the 
eigenvalue density of T2A |^, or by calculating 
the probability distribution function of (qq) [22\. 

Since the pioneering work of Nakamura 1 2J | , 
there have been several studies of TCQCD with 
/i 7^ using a variety of algorithms and actions 
I will briefly discuss some high- 
lights. Fig. SPshows strong couphng results 
for (xx): IT'S and {qq) vs. fi together with xPT 
predictions (^, showing acceptable agreement 
until saturation artifacts creep in for fi/rriTr ^ 0.6. 
Fig. ^ shows (xx)/(xx)o vs. rescaled chemical 
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Figure 3. (xx)j {qq) and vs. /i for /3 = 0, 
raa = 0.2, j/m — 0.1 on 4"* and 6^. 



potential for three different bare quark masses 
in TCQCD with adjoint quarks ||l|,|5|, together 
with the prediction (H) . In this case xPT appears 
to work well up to a; « 2 over a decade of quark 
mass. Since {qq) 7^ breaks a global symme- 
try, we anticipate a diquark Goldstone boson for 
fi > fj,o whose mass should vanish as j — > [Tl] ]. 
Physically this results in long-ranged interactions 
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between vortex excitations in a rotating super- 
fluid, and in propagating waves of temperature 
variation known as second sound. Simulations 
with j = slow dramatically in this regime due 
to the profusion of small eigenvalues of M [ p8[ . 
Fig. H shows the results of simulating the Pfafhan 
weight ( p^ with j '7^ on a much larger system 
|20[ , which enable the identification of a massless 
mode in the zero source limit. Finally, studies of 
the high density regime at T > have shown a 
strong first-order transition restoring the normal 
QGP state {qq) = (2§]. 

Although the TCQCD measure is obviously 
real, we have not discussed its positivity. In fact, 
it is possible to prove it positive for all cases ex- 
cept an odd number of adjoint staggered quarks. 
For A'' = 1 adjoint flavor, however, the sign prob- 
lem returns jl^. Additionally in this case, the 
Goldstone count does not include any baryons, 
and the superfluid condensate (qqs) in is for- 
bidden by the Pauli Exclusion Principle. Both 
facts invalidate the use of xPT. Indeed, the sim- 
plest local diquark operator that can be written 



is gauge- variant; 

qqlc = \ W^'t^ H- xi'x*'-] G 3 of SU(2), (12) 

leading to the possibility of color superconductiv- 
ity a la Georgi-Glashow in the high density phase. 
Simulations using the two-step multibosonic al- 
gorithm have been performed which reweight the 
ensemble taking the sign into account [|l^. Once 
/i > /io the sign starts to fluctuate, and by 
/i — 0.38 at TO = 0.1 (sgn(detM)) has fallen to 
0.30(4) on 43 X 8 m. Remarkably, once the ob- 
servables are also reweighted, all signals for the 
onset transition disappear, as seen in Fig.|. This 
is in accord with the expectation that there is now 
a separation between the lightest baryon mass 
and the Goldstone scale, and demonstrates that 
sgn(det) plays a decisive role in determining the 
ground state. Fig. ^ is probably the most expen- 
sive simulation of nothing happening to date. 

4. Flatland NJL 

The Nambu ~ Jona-Lasinio model has a long 
history as an effective description of the strong 
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interaction, and also has thermodynamic apph- 
cations . The Lagrangian is 

C = V^(^+^i7o+m)V^-^ [ii^^ f - (^75T^)'] (13) 

At zero chemical potential and for m = it has a 
SU(2)L(giSU(2)ij(X)U(l)B global symmetry, which 
for sufficiently strong coupling g'^ > sponta- 
neously breaks to SU(2)/(giU(l) b accompanied by 
the dynamical generation of a constituent mass 
E = g'^{ip'ip). Our current interest is the model 
in 2+1 dimensions, since apart from the obvi- 
ous computational saving there is an interact- 
ing continuum limit at g^ — > g"^, T,a p8| . 
For /i ^ 0, there is a strong first-order chiral 
symmetry-restoring transition at /ic — E j29| , 
which is completely separate from the Goldstone 



' , ■ appears to coincide ||31| with 
an onset transition separating the vacuum from 
a regime with ns o: /j,^, suggestive of a two- 
dimensional Fermi surface with Ep ~ /i in the 
chirally-restored phase. This is natural since the 
lightest baryons in the model are fermions. 

An obvious question is whether U(1)b is spon- 
taneously broken for ^ > fic by a superfluid 



condensate {x*^T2x) (note that we have 
surreptitiously slipped back into the notation of 
staggered fermions; the continuum translation is 
given in |^^). Superfluid condensation in this 
model would occur via a BCS instability at the 
Fermi surface, as in ^He. To investigate this we 
have performed Pfaffian simulations with diquark 
source term 1331 



j±qq± = j±(x*''t-2X ± 



(14) 



measuring both condensate {qq+) and associated 
susceptibilities x± = J2xil1±(^)l1±i^)) ■ Analo- 
gous to the axial Ward identity, we have 



X-(J- =0) = 



(15) 



The results are unexpected: fig. |^ shows a log- 
log plot of {qq+) vs. j+ from data from volumes 
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up to 32"^ extrapolated to the zero temperature 
limit Lt "> oo. The high density data {fXc — 0.65 
for our choice of coupling) suggest a power-law 
equation of state 



{qq+) ccf 



(16) 



with a = a{p,) falling in the range 0.2 - 0.3 for the 
values of fi examined. This is reinforced by the 
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susceptibility ratio R = \x+/x- \ plotted in fig. H 
using (p^ ) and (|l^) it is easy to show that R{i) 
should take the constant value a. The plateaux in 
fig. 1^ display this behaviour, with approximately 
the same a(/x) as those from direct fits to (|l6|). 
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There is a beautiful topological argument |35| 
that the only way to change the circulation k = 
§ Js-dl around a periodic volume is to create a 
vortex - anti-vortex pair and translate one of 
them around the universe in the perpendicular di- 
rection until they reannihilate, thereby increasing 
K by a quantum 2nKs/L. Since, however, the en- 
ergy required increases logarithmically with pair 
separation, the circulation is metastable, thereby 
demonstrating superfluidity. 

Why, then, does the NJL exponent 5{^) « 3 — 5 
differ from the XY value d{T) > 15? Stan- 
dard dimensional reduction does not apply since 
the limit Li oo is needed, implying that the 
fermion modes do not decouple. Further insight 
is gained from studying the fermion spectrum 
via the Gor'kov propagator Q — |Q. To 

probe the Fermi surface, we need to analyse non- 
zero momenta k; indeed the mass gap extracted 
from the decay of G{k) in Euclidean time ini- 
tially decreases as k y kp, corresponding to va- 
cating holes in the Fermi Sea, before rising for 
k > kp indicating particle excitations. The re- 
sulting quasiparticle dispersion relation E{k) is 
shown in Fig. ||. The detailed form of E{k) in- 



The strange behaviour of figs. is strongly 
reminiscent of the low temperature 2d XY model, 
which exists in a critical state for a range of T 
with continuously varying exponents 5{T), ri{T), 
which we can define analogously (the arrow de- 
notes a 2d vector): 

{qq)(xp ; (gq(O)gq(x)) oc (17) 

We therefore conjecture that NJL2+1 in its large- 
jjL phase describes a 2d critical system. The con- 
densate is washed out by long-wavelength fluctu- 
ations as J — > 0, but long range phase coherence is 
maintained via (p7|). Intriguingly, the most pre- 
cise tests of the 2d XY universality class are from 
experiments performed on thin films of superfluid 
''^He iQ. The superfluid current is related to the 
phase 6{x) of qq{x) via 
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Figure 9. E{k) for both free and interacting 
quasiparticles at /i = 0.8. 



(18) 



dicates a value for kp ^ and a Fermi velocity 
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vp — dE/dk\k=kF — 0-7c. This is characteris- 
tic of a relativistic Fermi hquid with a repulsive 
interaction between quasiparticles with parallel 
momenta p^ j . Most importantly, there is no evi- 
dence for a BCS gap A 7^ 0; certainly A ~ S can 
be excluded. The origin of the new universality 
class may therefore be attributed to the presence 
of massless fermions at the Fermi surface. The 
overall conclusion is that NJL2+1 describes a rel- 
ativistic thin film gapless superfluid. 

5. A conjecture about superconductivity 

All the systems we have been able to study at 
high density (ie. those with a real measure) ex- 
hibit superfluidity. What are the prospects for 
studying the breaking of a local symmetry? We 
have already seen that the only variant of TC- 
QCD with a potentially superconducting solu- 
tion, namely = 1 adjoint staggered flavor, is 
afflicted with a sign problem. The most promising 
microscopic model of high- Tc superconductors, 
namely the Hubbard model with on-site repulsion 
away from half-filling, also has an intractable sign 
problem in precisely the regime of interest ]37| . 
A more familiar example in a particle physics 
context is technicolor, in which condensation of 
quark pairs from different representations of a 
gauge group force its dynamical breakdown; since 
this requires chiral fermions, in complex represe- 
nations, a sign problem of some sort seems in- 
evitable. An exotic 2 + Id example is "ra-QED" 
p8[ |. This exhibits planar superconductivity by 
generating an electromagnetic photon mass by 
mixing with a second "statistical" photon via 
a mixed Chern-Simons interaction generated by 
fermion loops; for a non-zero CS coupling dynam- 
ical mass generation through a gauge-invariant 
condensate (iptp) 7^ is required. In Euclidean 
space, however, the CS term is pure imaginary 
and the resulting effective action complex; this 
also follows from the bare action since in this case 
the coupling iip^jsj^ip implies {75,^} 7^ 0, in 
turn implying detM 7^ detM*. Finally, of course, 
there is QCD itself. 

The following conjecture suggests itself: any 
system which exhibits the spontaneous breaking 
of a local symmetry by a pairing mechanism has 



a sign problem when formulated in terms of local 
gauge covariant degrees of freedom. 

6. Summary 

Significant progress - both technical and psy- 
chological - has been made in QCD(/i), result- 
ing in the first non-trivial LGT prediction in the 
(/i, T) plane. For once we have been lucky; the 
high-T low-// region where simulations can probe 
is precisely the regime of direct relevance to RHIC 
phenomenology. I anticipate much activity in this 
area in the coming year. At T = a quanti- 
tative description of nuclear or quark matter re- 
grettably seems as elusive as ever - however there 
are at least two model systems where LGT sim- 
ulations can now be said to be doing condensed 
matter physics ab initio, ie. with matter formed 
from the fundamental quanta of the theory. Two 
Color QCD, a confining theory with no Fermi sur- 
face, describes the condensation of tightly-bound 
bosons, and thus resembles superfiuid "^He. The 
NJL model, a theory without confinement but 
possessing a Fermi surface, displays unexpectedly 
interesting behaviour in 2 -I- Irf and certainly has 
the potential to exhibit a fully-fledged BCS mech- 
anism in 3 -|- Id, providing a relativistic analogue 
of superfluid '^He. Simulations of NJL3+1 will 
furnish non-trivial tests of model calculations of 
color superconductivity . Ultimately, though, 
a microscopic description of the superconducting 
state may actually require a sign problem. 
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